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A set of mode-coupled excitation equations for the slowly-growing amplitudes of dielectric waveg- 

^ ' uide eigenmodes is derived as a description of the electromagnetic signal field of a high-gain free- 

' electron laser, or FEL, including the effects of longitudinal space-charge. This approach of describing 

^\ ' the field basis set has notable advantages for FEL analysis in providing an efficient characterization 

' of eigenmodes, and in allowing a clear connection to free-space propagation of the input (seeding) 

' ^ ' • and output radiation. The formulation describes the entire evolution of the radiation wave through 

_ ^ ' the linear gain regime, prior to the onset of saturation, with arbitrary initial conditions. By virtue 

^ of the flexibility in the expansion basis, this technique can be used to find the direct coupling and 

^ ■ amplification of a particular mode. A simple transformation converts the derived coupled differen- 

■ tial excitation equations into a set of coupled algebraic equations and yields a matrix determinant 
• 1-H ' 

C/3 . equation for the FEL eigenmodes. A quadratic index medium is used as a model dielectric waveg- 

^ ■ uide to obtain an expression for the predicted spot size of the dominant system eigenmode, in the 

Oh< approximation that it is a single gaussian mode. 

ff-J . PACS numbers: 41.60.Cr, 42.25.Dd, 42.60.Jf, 42.50.Tx 



I. INTRODUCTION 



, The optical guiding of light in free-electron lasers (FELs) is a well-known phenomena that results during anipli- 

00 , fication when the coherent interaction between the source electron beam (e-beam) and the electromagnetic (EM) 
O • . 

^ . field introduces an inward curvature in the phase front of the light, refracting it back towards the lasing core of the 

^ ' e-beam . During the gain process the e-beam can behave like a guiding structure that suppresses diffraction, 

^>S ■ . 

^ reducing transverse power losses and enhancing the EM field amplification. For a sufficiently long interaction length, 

a 

■ - - ' the guided EM field eventually dynamically settles into a propagating, self-similar eigenmode of the FEL system, 
which, although describing a system in which the field amplitude grows exponentially as a function of longitudinal 
coordinate z, has a fixed form for the transverse profile and stationary spot sizey,(5|. 

Several different approaches currently exist for describing the FEL linear gain process. Analytic derivations of 
guided FEL eigenmodes have beenperformed previously by direct derivation of the eigenmode equations from the 
coupled Maxwell- Vlasov equationsP, Q; through solutions to the Maxwell equations using Green's functions Q|; or 
through expansion of the FEL signal field in terms of eigenmodes of a hollow, conducting-boundary waveguideQ], 
eigenmodes of a step-index fiber ^Sj], and free-space paraxial waves Such analytic descriptions of the 3-D FEL 
equations, particularly those that exploit radiation mode expansions to find the FEL eigenmodes, have certain utility 
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for providing physical insight into the character of the radiation. Since, in a high-gain FEL, the e-beam that propagates 
through the periodic undulator operates simultaneously as an optical source and as a waveguide structure, an expansion 
mode description of the FEL light permits investigation of the amplification, guiding and coupling efficiency of the 
individual basis modes to the e-beam. 

In general, the choice of analytical model and expansion basis depends on the FEL geometry and operational 
characteristics. Since the fields inside the e-beam are optically guided during high-gain, a guided expansion basis is a 
natural choice, but may be plagued by inefficiency in describing the EM field structure if the guiding characteristics 
of the e-beam are markedly different from those of the virtual waveguide that yields the particular basis set. (The 
waveguide is referred to as 'virtual' because no such waveguide is assumed to be present in the physical system.) In 
Ref [5| an analysis is presented for calculating the gain-guided FEL supermode by means of a field expansion using 
eigenmodes of a hollow waveguide with conducting walls. This method can be useful at short wavelengths (e.g. x-ray 
FEL), even when no such physical boundary exists and the optical beam is guided only by the e-beam. Clearly, 
in this example analysis and in general, this approach is legitimate if the presence of the artificial boundary does 
not affect the physical result. The virtual waveguide dimensions must be taken wide enough so that the fields in 
the actual problem vanish at the virtual boundaries, but not so wide that many modes are required to describe the 
radiation field such that the calculation is inefficient or does not converge. As an alternate approach described here, 
a general method that uses eigenmodes of a smooth, transversely anisotropic dielectric waveguide can be used to 
model the field while avoiding both the undesirable infiuence of artificial boundary conditions and precise a priori 
knowledge of the characteristic transverse dimensions of the waveguide features of the e-beam. This approach also 
provides flexibility in the form of the expansion basis, which is determined by the refractive index distribution. Of 
particular interest in any expansion mode approach is the coupling to the naturally occurring Hermite-Gaussian (HG) 
or Laguerre-Gaussian (LG) modes that describe free-space wave propagation in the paraxial limit. This connection 
is useful for relating the FEL modes to free-space propagating modes, which are present both before and after the 
exponential gain process occurs in the undulator. With this connection one may, for example, compactly and clearly 
describe the input radiaiton needed for efficiently seeding the FEL. One may also robustly describe the propagation 
of the FEL light after saturation, allowing a clear understanding of the demands made on downstream optics, and on 
expected mode structure at FEL diagnostics or experiments themselves. In the virtual dielectric waveguide expansion, 
a form of the refractive index that varies quadratically in the transverse dimension yields a basis of guided LG or 
HG paraxial waves. This facilitates the desired natural description of the FEL radiation, and is the motivation of the 
present work. 

Optical guiding in an FEL occurs in the exponential gain regime, when the focusing effects in the source e-beam tend 
to balance the natural diffraction of the generated EM radiation. Obviously, an expansion set consisting of guided 
modes can be used to accurately describe the field during long sections of gain-guiding, over which the otherwise 
free-space fields have many Rayleigh lengths to diffract. Such an expansion can also accurately describe short sections 
of dynamic profile evolution (like the FEL start-up period in a high-gain amplifier). But guided modes may not 
efficiently capture (in the sense of providing a compact description) the field behavior during long sections of weak 
guiding when diffraction is dominant, like during saturation or low-gain. The regimes of validity for a description 
using guided modes may be estimated by inspection of the relative values of the optical Rayleigh length Zfi = kr^ 
(assuming for the moment that the field has the same characteristic radius as the e-beam, rg), the e- folding radiation 
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power gain length Lq, and the overall length of the FEL interaction length L. In the guided mode expansion here, 
we restrict our attention to high-gain {Lq ^ L) FEL systems, prior to the onset of saturation, for which the weakly 
diffracting [Lq <^ L < zji) or strong-guiding [Lq < zj^ < L) conditions are valid, and the guided mode expansion 
description can be used efficiently. 

In this work, a virtual dielectric expansion description of high-gain FEL radiation is presented. The excitation 
of the slowly-growing mode amplitude coefficients in the presence of a local source current is derived. The e-beam 
is linearized in the cold-beam fluid approximation and a set of coupled excitation equations, modified from Ref 
is obtained. The coupled equations are then recast by a simple transformation into a matrix equation for solutions 
of the FEL supermode as a function of parameters for energy detuning, longitudinal space-charge effects and mode 
coupling to the e-beam. The theory presented here provides a framework for detailed numerical studies performed 
in a follow-up paper (II. Modeling and numerical simulations). There, results that highlight the advantages of the 
dielectric eigenmode expansion approach are presented and analyzed. 



II. DIELECTRIC WAVEGUIDE EXPANSION MODES IN THE PRESENCE OF A LOCAL SOURCE 



The radiation fields ^(r, t) and Hjr, t) emitted by the currents in the FEL can be expanded in terms of transverse 
radiation modes {£_j^,2ii_} of ^ guiding structure, with slowly-growing amplitudes that vary only as a function of the 
axis and e-beam propagation coordinate, z. Neglecting backward propagating waves (which is well justified in an 
FEL), the general field expansion is, 



(1) 



where 2i_Lg — {^/Zq)ez x is the magnetic field expansion mode, kzq is the q*"^ mode axial wavenumber and the 
impedance is Zq — {k/kzq)\/ fJ^o/^o for TE modes. In general, the sum extends over both the guided and the cutoff 
modes, and the modes form a complete set[lo|. The expansion modes are orthogonal and normalized to the mode 
power 



VqSq,q' = ^Rc 



[l^qirAxnlq'ir±)]-ezd'r^ 



(2) 



In free-space, the equation for the full time-harmonic fields in the presence of a source is 



(3) 



where the total field is E_{r,t) =Re[£^(r)e^*'^*] and k — uj/c is the free-space wavenumber. The transverse charge 
density gradient of the e-beam in an FEL is typically small and is neglected. For a given choice of expansion basis, the 
field expansions in Eqs. ([!]) can be inserted into Eq. ([3]) to obtain a differential equation for the amplitude coefficents 
(^9(^)151. Here, we explore the case where the expansion mode £±q is an eigenmode of a dielectric waveguide with 
refractive index n{r±). Assuming small transverse variation, 'Vn{r±)'^ <C k, the eigenmodes can be regarded as 
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dominantly transverse, and the dielectric eigenmode equation is 

Vi£^,(rx) + Hr^fe - kl^]£^^{r^) - 0. (4) 

The use of dielectric eigenmodes as an expansion basis must be accompanied by additional terms in Eq. ([3]) . These 
terms are included to offset the virtual polarization currents that arise from the refractive index of the virtual 
waveguide, since no such waveguide exists along the e-beam axis in the physical system. The polarization current is 
given by J_p{r) — iujeox(r±)E^{r±) = iijjeo{n{r±)'^ — ^)K{t±) where n{r±) = \/e{r±)/eQ. Mathematically, this process 
amounts to adding (n(rj^)^ — l)k'^E_{r±) to both sides of Eq. ([3]). From the dielectric eigenmode equation the 
modes are assumed to be dominantly transverse, and Eqn ([3]) becomes 

y^E^ + n{ri_fk'^E^ = -iLu^ioli_ + - l)fc^l(ri). (5) 

Plugging in the expansion fields from Eq. ([T]), the dielectric eigenmode equation in Eq. (j4]) is used to eliminate the 
transverse Laplacian term, and the excitation equation for the mode q in the presence of a local source current is 
given by, 

where 

<y ^wjl t''^''^^' ~ • lUr±)d'r^- (7) 

The difference between the axial wavenumbers of the modes is Ak^qq' = k^q — k^q' . The term ^, characterizes 
the mode overlap in the dielectric, and physically represents the virtual polarization currents that are necessarily 
subtracted when using eigenmodes of a virtual dielectric waveguide. 

III. ELECTRON-BEAM FLUID AND COUPLED EXCITATION EQUATIONS 

A linear plasma fluid model for a cold e-beam (negligible energy spread) can be used to describe the signal excitation 
in an EEL interaction^. A relativistic e-beam in an EEL experiences transverse oscillations, or wiggling, driven by 
an interaction with a periodic structure or undulator. This periodic motion drives an axial ponderomotive force 
that modulates the axial electron velocity such that, to first-order, the axial velocity of a cold beam within a static 
undulator can be expanded as Vz{r,t) — Vzo + Re['D2i (r)e^*'^*] where Vz„ — PzC is the d.c. component and Vz^ is the 
modulation at signal frequency lo. Longitudinal variations in the velocity like those found in planar undulator systems 
are ignored for the moment. The velocity modulation Vz-^ develops a density bunching modulation that is similarly 
described in a linear model as n(r, t) — no/(?'j^) -|- Re[ni(r)e~'"*] where Pq is the on-axis electron density and f{r±) 
is the transverse density profile of the e-beam. The lowest order a.c. component of the longitudinal current density 
results from both the axial velocity and density modulations and is Jz{l) — ~e[no/(r'_L)i)zi (n) + Vz„f\i{r)]. If the 
transverse divergence of the current density modulation is assumed small Vx • Ji_ ^ dJz/dz, the continuity equation 



5 



can be written as 

d 



, Jz = -«weni(r) (8) 
az 

The relativistic force equation for the axial vefocity modulation is, 

^~ (^) _ ^ ^ ^ ^ 5)^] (9) 

where the axial space-charge field E^'^ — Jz/iujeo is due to the current density modulation. It is assumed that the 
e-beam radius is large compared to the bunching wavelength in the e-beam frame, tq > X^z, so the fringing fields 
and the transverse space-charge effects are neglected. The interaction between the transverse electron motion and 
the transverse magnetic fields results in an axial ponderomotive field {v x B_)z = C!q{z)(S'pm,q{r±)e^'-^^'''^'^^^^ where 
iS'pm,q{r±) = ^[^±q ^ S.±tu + ^Itu ^ IS.±q] ' 62, B±w is the trausvcrsc magnetic field of the undulator and v±q{r±) is the 
transverse electron velocity due to the Lorentz force of the q*^ mode of the signal field. The undulator wavenumber is 
identified as — 27r/Au;, and the transverse velocity due to the magnetic undulator field is Vj_^ = {—icK/"/)ez x 
where K — e\B _\_w\ / mckw is the undulator parameter, and Cw = B_^^/\B±w\ is the unit vector of the undulator field. 

By combining Eqs. ([8]) and ([9]) and using the definition of the modulated longitudinal current density, the density 
bunching evolution can be expressed as a second order differential equation. We obtain the result from Ref. |5|] with 
transverse fields for the density modulation evolution during the FEL interaction: 



dz^ 2*-- 



Ml) = ^ \ - T^ikzq + fc^)C,(z)4,™,,(r^)e*('=-+'=™)^ (10) 



where the longitudinal relativistic plasma frequency on axis is ujp^ = e'^Uo/jj'^eQm. The beam current is Iq = 
— ewzo"o JJ f{T±)d'^^± and the effective beam area through the normalization condition JJ f{r±)d'^v± — A^. For 
a uniform cross-sectional distribution of the e-beam, ](r\_) = 1 for r < vq and zero otherwise. For a gaussian 
distribution, f{r±) =exp(— r^/rg). 

The transverse component of the current density that excites the signal wave is written in terms of the density 
modulation as 

l±{r) ^ -len,{r)v^^e-'''-\ (11) 

The charge density modulation ni(r) appears both in the field mode excitation equation ^ and in the density 
modulation evolution equation (jlOp . These equations are coupled through the ponderomotive field, which illustrates 
the relationship between the density modulation and the excited signal field in an FEL. Both equations can be 
simplified by expressing the density modulation as a sum over the expansion basis functions and slowly-varying 
amplitudes. 

Ml) = ^5]i?,(z)f\,(r^)e'^^ (12) 

9 
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Plugging this expression into Eqs. ((TT|) . ((TO)) . and and integrating over the transverse dimensions in Eq. (|10p . 
the orthogonahty of the eigenfunctions can be used to simphfy the equations and to write both coupled equations in 
terms of the slowly growing amplitudes Cq{z) and Bq{z) of the signal field and the density modulation, respectively. 
This yields the coupled FEL excitation equations: 



q' 



(13) 



where 9q — uj/vzg — {kzq + k^) is the characteristic detuning parameter for a given mode q, ^q — Kk"^ / A^yk^q, and 
9p ~ LUpg/vzg is the longitudinal plasma wavenumber of a uniformly distributed electron beam profile used in a ID 
model. The factor g* = (cz x e^i) • e* is the polarization alignment factor and measures the relative alignment of 
the transverse electron motion in the undulator (e^ x e^) with the (complex-conjugated) E-field mode polarization 
direction (e*). When the wiggling motion direction matches the mode polarization (which may not be the case for 
the input mode in a seeded FEL scenario), gq = 1. 

The coupling between the e-beam and the signal field is given by the mode coupling coefficient: 



Qq.q' = JJ^(^-«' + // /(^^)4>™.9'KHi^„ • £lqiri_)d:'r^ (14) 

where JJ= [Jo(a) — Ji{<^)]'^ is now included for a strong planar undulator (JJ=1 for a helical undulator geometry). 
Jq and Ji are the first and second order Bessel functions and a — K'^/{A + 2K'^). The ponderomotive field S^m.q' is 
evaluated explicitly for TE modes in the Appendix, and we obtain a simplified form, 

Q^ ^, = jj 02 (fc,2_±M!(^)\*^^,F,,,,. (15) 

where Fg.q/ is referred to as the beam profile overlap coefficient and quantifies the spatial overlap of the e-beam profile 
with the expansion modes q, g': 

JJ |f±g(»-±)Pd2rj. 

This coefficient also appears in connection with the longitudinal plasma wave dynamics (left hand side of the second 
equation in Eqs. I13p . 

It also may be illuminating to also define the current bunching amplitude: 

^«(^) = -T^e"^^"^""^''' / / L^_{'d-~Cq(r^_)d\^ (17) 



4-p, 
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which is interpreted as the slowly- varying amplitude of the transverse optical current. Replacing with it's definition 
from Eq. (jlip and plugging in the density modulation expansion from Eq. (|12p . the integration over the transverse 
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space is straight-forward by the orthogonahty of the eigenmodes. Thus, a simphfied expression is obtained for Xq{z) 
in terms of the density modulation amphtudes: 

= -i^,9;B,{z). (18) 

Substitution of this into Eq. (fT3|) allows one to solve the coupled FEL excitation equations in terms of the optical 
current amplitudes. 

In the ID limit, the axial wavenumbers are degenerate k^q ~ k, the e-beam profile is constant, f{rj_) 1, and 
the overlap factor becomes IFq.<j' = 5q,qi . The e-beam coupling coefficient is then written in ID simplified form, 
Q = 33k{9pK/2^/2'yPz)^ = {2kwp)^, where p = (JJe^ii'^no/32eo7^mc^fcJ)^/^ is the well-known Pierce parameter 
often used in FEL theory. 

The first of the coupled equations in Eqs. (fTS]) describes the excitation of the mode amplitude Cq of a dielectric 
waveguide eigenmode due to the transverse wiggling motion of the bunching current throughout the FEL interaction. 
The second equation in (|13p describes the evolution of the density bunching amplitude, which is excited in the e-beam 
by the EM signal field. The effect of the longitudinal space-charge in the beam is calculated in the second term of the 
second equation, and takes into account the effects of longitudinal plasma oscillations (Langmuir waves) in the FEL 
interaction. The effects of fringing fields due to a transverse variation in the axial space-charge field can be examined 
in a more complete 3D scenario ll', but are presently neglected in the approximation that rg > X'jz- 

The initial conditions for Eqs. (|13p specify the operating characteristic of the FEL. For example, when operating as a 
single-pass amplifier (seeded FEL) there is negligible initial density and velocity modulation {Bq{0), dBq(z) /dz\z=Q — 
0) and the initial seed field is non-zero (Cg(0) 0). For a self-amplified spontaneous emission (SASE) FEL, the 
amplified shot noise can be related to the pre-bunching conditions {Bq(0) ^ 0, dBq(z)/dz\z=Q ^ 0) and the input 
signal field vanishes (Cq(0) — 0). 

The coupled expressions in Eqs. p3|) describe the evolution of the e-beam modulation currents and the excitation 
of the EM signal field during the FEL interaction, inclusive of longitudinal space-charge effects. During high-gain 
operations when the radiation field grows exponentially, the self-guiding effects become dominant over the diffraction, 
and the radiation field can be accurately described by a collection of waveguide modes. The virtual dielectric eigenmode 
basis used for the radiation field expansion can be any complete set that satisfies the dielectric waveguide equation 
(l4|). In the general expansion presented here, the form for n{r±) can be freely chosen to yield different functional 
forms for the expansion modes, allowing fiexibility to choose a basis that is perhaps better suited to describe a given 
FEL system. Some of the most typical forms include an index with a quadratic spatial dependence which results 
in Hcrmitc-Gaussian or Laguerre-Gaussian expansion modes 12| (examined in detail for LG modes in part II of this 



work), or a step-profile optical fiber [3|, [10| which yields Bessel functions. The closer that the choice of the virtual 
dielectric waveguide distribution is to the real solution to the given FEL system, the fewer the number of modes that 
will be required to converge the equations in Eqs. (fT3|) to the correct solution for the FEL. 
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IV. SUPERMODE MATRIX SOLUTIONS 

In the high-gain regime of an FEL with finite e-beam width, the radiation field tends to concentrate near the beam 
and produce a power amphfied radiation wave with a self-similar transverse field distribution that propagates along the 
interaction length. This specific complex-valued combination of the expansion modes is referred to as the supermode, 
or the eigenmode of the high-gain FEL. In a long undulator, the fundamental supermode evolves spontaneously since 
it has the highest gain. The initial conditions (i.e. the transverse profile and phase distribution of the injected 
radiation field or density pre-bunching) affect the mode excitation composition and evolution. However, these initial 
conditions only affect the supermode establishment length. Evcnutally, in an undepleted system, the supermode will 
become dominant in a long enough undulator, prior to the onset of saturation. 

To find the optimal injection parameters that match and expedite the establishment of the supermode, one must 
solve Eqs. (|13p with the appropriate initial conditions. To simply find the system supermodes however, it is enough 
to find the eigensolutions to Eqs. (fT3|) . These are the combinations of the expansion mode profiles that propagate 
self-similarly, i.e., with constant amplitude coefficients and with distinct complex wavenumbersQ. In the presence 
of gain, each supermode wavenumber will be different from the wavenumber of free-space and can be written with a 
perturbative term dk that is due to the FEL interaction: 

ksM = k + Sk, (19) 

where Re{5fc} anticipates an effective modified refractive index to that of free-space, and Im{5fc} is related to the 
gain. Since the supermodes evolve after the initial startup period and have fixed transverse profiles along z, one can 
substitute Cq{z) — hqe^'^''^'"'^^"''^^ for the mode amplitude coefficients Eq. ([T]). The mode amplitude coefficients bq are 
constants, and the z-dependence is contained solely in the mode-independent exponential term. The time-harmonic 
electric field of the supermode is then 

EsM{r)=[Y,hqlq{r^)y'^--^. (20) 

9 

This describes a transverse field that is fixed in transverse profile, but is growing exponentially in amplitude along 
z. Inserting this transformation into Eq. (|13p converts the coupled second-order differential equations into a set of 
coupled algebraic equations: 

{5k - 9f [{5k - Akq)bq + J2 V] - E I l^^^'J- [(-^"^ - + E ] + E Qi^i'"^'" = 0' (21) 

g' j K ^'^ q" J q'" 

where 9 = to/vzf, — k — k^ is the detuning parameter for a ID model {kzq — k), and Akq = k^q — k. The coupled 
equations in (|2f p can be solved to yield values for the supermode coefficients hq in terms of the perturbation 5k. It is 
convenient to write Eqs (|2ip in a simplified matrix determinant form 
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The matrix elements of M are given by M, 



[kzq' /k^q)¥q^q>, and similarly for ^ = {^1^^,/}, Q_ = {Qq,q'}, and 



Afc = {/S.kq5q^qi}. The matrix I is the identity. 

The solutions to Eq. yield 3iV solutions for 6k, where N is the number of expansion modes. Each 5k can then 
be inserted in Eq. ((2T|) to find a non-trivial solution (if one exists) for the mode amplitude coefficients of an eigenmode 
of the FEL system. From Eq. ([20]) and the definition of ksM in Eq. (fT9|) . it can be seen that the solution for 5k with 
the most negative imaginary component drives the highest gain, and dominates over the rest of the eigenmodes. This 
value is used in solving Eq. (|2ip and will yield the coefficients of the dominant supermode, with the 3D power gain 
length, or e-folding length, given by Lq — (2|Im{5fc}|)~^. 

We note that when Afc, — 0, the matrix equation in ((22)) reduces to a generalized matrix form of the canonical 
FEL cubic equation with longitudinal space-charge. In the additional limit of a large, effectively constant transverse 
beam profile, f{r±) ^ 1 so ¥q.q' = 5q,qi and Qq,qi Q5qy, and Eq. (|22p takes the form of the familiar ID FEL 
cubic equation with the mode-independent beam coupling parameter Q and the axial plasma wavenumber 9p. 



During high-gain, the proper balance between the natural diffraction of the coherent radiation and the guided 
focusing of the radiation due to the e-beam determines the eventual spot size wsm of the EM supermode field. This 
can be obtained with the dielectric expansion formulation in a natural way, through the solutions to the excitation 
equations (jl3p or the supermode determinant equation (|22p . using N modes in the field expansion. Clearly, choosing 
a suitable expansion basis (i.e., one that is close in form to that of the supermode field) can greatly reduce the number 
of expansion modes required to precisely describe the FEL supermode. In some cases, it may even be sufficient to 
use only a single mode to model the field, in order to estimate the relevant supermode characteristics without solving 
the full equations. A single mode approximation can also be used to streamline the computation of the full solutions 
by providing an approximate value for the expansion basis waist size wq which is used as a scaling parameter in the 
basis eigenfunctions £q{r±;'Wo). The supermode waist size (and thus the parameter wq) can sometimes be roughly 
approximated as equal to the e-beam radius wsm ^ ^0 for an axisymmetric e-beam distribution, but the actual value 
may vary significantly depending on the operating parameters of the FEL. To better quantify this estimate for the 
general case, we analyze the guided supermode spot size by approximating it as a single gaussian mode (SGM) that 
propagates with a fixed spot size, Wg. Gaussian modes, like the HG and LG modes discussed previously, are found 
as solutions to Eq. ([4]) with a refractive index n(r)^ that varies quadratically in the transverse dimension, describing 
what is known as a weakly-guiding quadratic index medium, or QIM. In this context, the SGM analysis also provides 
insight into the physical origin of the characteristic guiding of the supermode in terms of guided paraxial gaussian 
modes from conventional fiber waveguide theory. 

The refractive index that yields axisymmetric guided gaussian eigenmodes in Eq. ^ is 



V. SINGLE GUIDED GAUSSIAN MODE APPROXIMATION 




(23) 



where zn — kwg 



/2 is the Rayleigh length and Wg is the characteristic rms radius of the gaussian field profile given by 
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£q{r±; Wg). The axial wavenumber, in the paraxial approximation is therefore given by 



2 

^9 



kM^kno — . (24) 



A single gaussian mode then corresponds to only the (q, q') = (0, 0) matrix elements in the supermode determinant 
in Eq. ([22|). For np ~ 1, an explicit analytic expression is obtained for the dielectric mode coupling parameter 
Kq Q = — 1/fcw^. Ignoring the effects of space-charge waves {9pL < tt), the form of Eq. for a single gaussian mode 
is 

{Sk - [5k + + Qo.o - 0. (25) 

where Afco — —2/kw'^ has been obtained from Eq. (|24p . Equation (P5)l then becomes a modified ID EEL cubic 
equation where Qo.o = Q ^o,o is the SGM gain parameter, which describes the modification of the ID e-beam mode 
coupling parameter Q by the filling factor (Eq. [TB)) of the fundamental gaussian EM mode, Fo,o- One can then define 
a shifted perturbation parameter: 



Sk = Sk + -^ (26) 
kw^ 
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and then set the detuning to 6 = —1/kWg. Equation then becomes a simple ID FEL cubic equation at resonance, 

Sk^ + Qo,Q = (27) 

In this form, the solution for 6k that corresponds to the dominant high-gain mode is straightforward, and well-known 
from FEL theory: 

5k = ^-^Ql'l (28) 

In a high-gain FEL, the supermode wavenumber ksM in Eq. (jl9p is defined as the wavenumber of a plane wave k 
that is modified through the FEL interaction by the perturbative factor 5k. Combining the definition of fc^^v/ in Eq. 
(fro)) with 5k from Eq. we can write 



ksM ^k + 5k--^ (29) 
kw^ 



a 



for a single gaussian supermode. Since the supermode is assumed to propagate as a guided gaussian mode with a 
fixed spot size and velocity, the real (propagating) part of the supermode wavenumber is equated to the propagating 
gaussian mode wavenumber in the dielectric waveguide: 

Re{fcsAf} = fc.o. (30) 
Combining the real parts of Eqs. ([28| . (|29|) and ([30]) . and recognizing no = 1 + Q^/^ /2k ~ 1 as the effective refractive 
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index on axis from a ID FEL model, we obtain a relation for the supermode spot size in terms of the ID coupling 
and the filling factor 

i^ = V<>-C)- (31) 

Assuming that the e-beam has a gaussian transverse profile f{r±) =exp(— r^/rp), the filling factor is simply Fq^q = 
{l + 'Wg/2rQ)~^ . We finally obtain an expression for the supermode spot size Wg in terms of the ID coupling parameter 
Q and the e-beam radius rot 

r 2 1 3 

(1 + ^)1-77:27^ =1 (32) 



'0 



Solutions to Eq. (j32p for Wg can be easily obtained numerically for given FEL running parameters. The calculated 
SGM spot size provides a useful prediction of the supermode spot size wsm for a gaussian e-beam lasing at the 
fundamental and which in turn suggests an appropriate value for the expansion waist wq to facilitate computational 
efficiency for solving the full excitation equations. 

It is useful to note that while ksM was defined as a modification of a plane wave in Eq. (fT9|) . in the SGM 
approximation, one could also consider defining ksM through a modification of a free-space paraxial gaussian mode. 
In other words, to consider the effect of the FEL interaction on converting a free-space gaussian wave to a guided 
gaussian mode of the FEL. 

Free-space modes evaluated at the optical beam waist wq have an associated axial wavenumber 

^zo = — (33) 
Expressing the supermode wavenumber as a modified paraxial wave of free-space, it takes the form 

ksM ^ kzo +Sk (34) 



where, just as before, the supermode is a guided eigenmode of the QIM, Iie{ksM} = kzo, but it is now also expressed 
as a deviation from the gaussian free-space mode kzo by the FEL interaction given by 5k. Equating the real parts 
from Eqs. (|19p and (p4|) . we find that k + Ke{Sk} = k^o + Ke{dk}. By rearranging and using k^o — fc ~ —1/kWg for 
a single free-space mode, we recover the exact expression for Sk in Eq. (P^ . Thus, the SGM approximation of an 
FEL with a gaussian transverse e-beam profile describes the conversion of plane wave into a guided gaussian mode at 
a detuning of = — 1/ kwg in an equivalent way as it describes the conversion of a free-space paraxial gaussian mode 
into a guided gaussian mode, operating at resonance {9 — 0). 



VI. CONCLUSIONS 



The signal field of a high-gain FEL has been described in terms of a sum over orthogonal eigenmodes of a virtual 
dielectric waveguide. A set of coupled excitation equations for the e-beam density modulations and the field amplitudes 
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in the presence of longitudinal space-charge effects has been derived. This approach can be effectively used to predict 
relevant FEL parameters such as signal field spot size and intensity distribution for any point along the interaction 
length, and for any arbitrary initial conditions. It further provides a novel and practical method to analyze the coupling 
and generation of particular modes, such as those associated with free-space propagation, in the FEL system. It is 
a general method for calculating the optically guided FEL eigenmodes for any transverse e-beam current and profile 
distribution, based on a straightforward numerical solution to an algebraic dispersion equation. A follow-up analysis 
of this work will explore these concepts in detail. 
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APPENDIX A: E-BEAM COUPLING COEFFICIENT 



A generalized, explicit expression for the mode coupling coefficient Qq^q' of any spatial expansion basis can be 
obtained starting from the ponderomotive field 

1 -* 

<^pm,q{ri_) = -[V-^q X B^ru+V-lw X ^±g] ' ^z- (Al) 



From the relativistic Lorentz force equation we can write the transverse velocities as 11 



(A2) 



where in the last equality for v_^^ it has been assumed that £_^^ = —{Lj/kzq)ez x ^A_q for TE modes, in keeping with 
the dominantly TEM mode expansion sets of the present work. With the equations in (|A2|) the expression for S'pm,q 
becomes 

The complex phasor notation for the undulator field =Re{&„,e^''^"'^} allows a general description for any type 
of undulator polarization. Defining a general transverse field polarization unit vector as — M.w/\^w\, for a linear 
undulator = iy (with vertical polarization) , the field is given as 

= \Bw\ey cos z ( A4) 
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and in a helical undulator = {e-x ± iey)/\/2 for right-handed (+) or left-handed (-) field orientations and 

= (e^; cos kwZ ± By sin z) . ( A5) 

v2 

The undulator parameter is defined as K = e\Bw\/{kw'mc). Note that the same K value corresponds to an on-axis 
maximum field amplitude of \Bw\ in a linear undulator, and an on-axis maximum field amplitude of |,Biu|/\/2 in a 
helical undulator. In terms of Eq. (jASp can be written as 



= ^^^^^(^^ X -^^.K). (A6) 

Substitution of this expression and the mode power normalization from Eq. ([2]) into Qq^q' results in a general expression 
for the mode coupling coefficents: 



(-) 9l9,'^,,,- (A7) 



where cjq = (e^ x e^) • is the polarization alignment factor, which measures the relative direction of transverse 
electron motion in the undulator compared with the electric field polarization of the mode q. The spatial overlap 



factor Vq,qi is defined as 



//|Ex,('-x)P<i^i-x 



The general expression for the mode coupling coefficient in Eq. (jA7P permits analysis of mode coupling for radiation 
modes that vary transversely in their complex field amplitude and in their polarization relative to the undulator. It is 
also clear that in order to obtain maximal coupling between two modes q and g', they must be polarization matched 
to each other (g^ — (jqi) and to the direction of electron motion in the undulator {gq — 1), as well as have finite spatial 
overlap with both the e-beam distribution. In the single mode limit [q ~ q'), Eq. (jASp simplifies to the commonly 
used "filling factor" parameter. 
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